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ABSTRACT 

The well known scaling laws relating critical exponents in a second order phase tran- 
sition have been generalized to the case of an arbitrarily higher order phase transition. 
In a higher order transition, such as one suggested for the superconducting transition in 
Bao.6Ko.4Bi03 and in F^S^CaC^Os, there are singularities in higher order derivatives 
of the free energy. A relation between exponents of different observables has been found, 
regardless of whether the exponents are classical (mean-field theory, no fluctuations, in- 
teger order of a transition) or not (fluctuation effects included) . We also comment on 
the phase transition in a thin film. 

1. Introduction 

In a recent letter (Kumar, Hall and Goodrich 1999; Hall et. al. 2000), there is a proposal 
that the superconducting phase transition in Bao.6Ko.4Bi03 (BKBO) is of order four, in the sense 
defined by Ehrenfest (1933). Thus the first three derivatives of the free energy with respect to the 
temperature T, the entropy S, the specific heat C and the temperature derivative of the specific 
heat dC/dT are continuous across the phase transition. As are the first three derivatives with 
respect to the magnetic field H; the magnetization M, the susceptibility x an d dx/dH. The fourth 
order derivatives, namely d 2 C/dT 2 and d 2 x/dH 2 are discontinuous. The discontinuity may even 
be a small power-law singularity such as a A-point anomaly but in d 2 C/dT 2 and d 2 x/dH 2 . 

There are several materials (Kumar, unpublished 2001) which are possible candidates for a 
higher order phase transition (see the discussion at the end of this paper). One class consists of 
Barium and Bismuth based cubic perovskites, doped with K, Rb or Sr. One of the earliest anoma- 
lies about the original material, Ba(Pbi_ x Bi x )03 (BPBO) was an absence (Methfessel, Stewart, 
Matthias and Patel 1980) of the usual specific heat discontinuity at the superconducting phase 
transition. Many of the materials listed in (Kumar, unpublished 2001) have not been studied for 
their specific heat yet, nor for any anomalous feature in the associated magnetization. For a second 
order phase transition, the phase boundary in the H-T phase diagram is described by, 
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where T c is the transition temperature. If AC is zero then either the phase boundary is flat, or Ax 
must also be zero and the transition is of a higher order. In most materials which are paramagnetic 1 , 
it is quite unlikely that Ax = 0. Most of the materials discussed above however are diamagnetic 
in the normal state. Here the condition Ax = may be more easily satisfied, thus making them 
possible candidate materials for a higher order phase transition. 

Following an extensive analysis of the temperature dependence of the specific heat near the 
superconducting transition, Junod, Erb and Renner (1999) have proposed that the superconducting 
transition in Bi2212 (E^S^CaC^Os) is similar to the transition one encounters in Bose-Einstein 
condensation in an ideal Bose gas. From the specific perspective here, this is an example of a third 
order phase transition. We elaborate on this topic below in the discussion section. 

The purpose of this paper is to further explore the thermodynamics of a p th order phase transi- 
tion. Let us define the terms: the condensed state free energy —F(T) ~ t p ~^, where t = (1 — T/T c ), 
p is an integer and /i / 0, a small correction to p. In general, fi should also have a subscript p 
signifying its dependence on p. We will often omit this subscript for convenience. The term classical 
will refer to the case ji = and the term non-classical will include effects of a non-zero fi. For 
example, for a second order transition p = 2 and \x = a, where a is the specific heat exponent. 
The term non-classical then encompasses all fluctuation effects related to the ground state. In the 
results reported in refs. 1 and 2, the exponent for the temperature dependence of the free energy is 
3.6. The suggestion thus being that p = 4 and fi = —0.4, the specific heat exponent a = —1.6 < 0, 
i.e., there is no divergence in the specific heat but that the second derivative of the specific heat 
d 2 C/dT 2 has a A-like divergence with an exponent fi. Similarly, there may be a A-like divergence 
in d 2 x/dH 2 with a corresponding exponent. 

We have focused the discussion here to a free energy depending on magnetic field and tem- 
perature. In general, there are other mechanical variables which can be accommodated by simply 
replacing the thermodynamically conjugate pair (M, H) by the appropriate combination (for ex- 
ample volume V and pressure P). Moreover, we will also discuss two features which are specific to 
a superconductor. In a superconductor, the magnetic field influences the charge motion via a gauge 
coupling. Indeed, if the Cooper pairs are assumed to be in a spin singlet state, as we assume in 
the following, then the magnetism is mostly due to the orbital contribution and the susceptibility 
has a special temperature dependence. The scaling exponents are then specific to that case. The 
second superconductivity feature here is the existence of two different critical fields. In the H-T 
plane, flux expulsion occurs at H c i(T), which is lower than the upper critical field H C 2(T) where 
superconductivity is destroyed. Both of these critical fields vanish at T c . In a second order phase 



1 In the Abrikosov state, the magnetization is negative but the susceptibility is paramagnetic and is usually larger 
than the normal state paramagnetism value. Thus one often finds a discontinuity in \ at H C 2- In contrast, the 
materials discussed here are often diamagnetic. It is easy then to smoothly connect the magnetization across the 
superconducting transition so that around the upper critical field H C 2, the susceptibility is zero on both sides of the 
transition. The discontinuity is then in the nonlinear susceptibility as expected for a higher order transition. 
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transition they are both linear in reduced temperature and their ratio, the Landau parameter k, is 
a constant, independent of temperature. In a higher order transition (Kumar, Hall and Goodrich 
1999) k is temperature dependent, diverging at T c . The lower critical field also is a measure of the 
superfluid density through the London penetration depth. 

Thus the well-known scaling laws appropriate for a second order phase transition can be 
generalized for an arbitrary order phase transition as described in Sec. II below. In Sec. Ill, 
we study a relationship between the specific heat and the London penetration depth, which is a 
relation between the exponents p — \x and an exponent for the temperature dependence of the lower 
critical field H C \{T). Section IV contains a discussion of finite size effects. In this section we also 
discuss a Kosterlitz-Thouless type phase transition, the binding-unbinding of a vortex-antivortex 
pair, as well as the irreversibility effects associated with the melting of a flux lattice. Finally, Sec. 
V contains a summary of our results including a discussion of possible candidate materials and field 
theory models for a higher order phase transition. 



2. Scaling Laws 

Scaling laws are consistency checks based on the (magnetic) field-temperature dependence of 
the experimental observables. The free energy is considered as a function of temperature T and a 
mechanical variable, say, the magnetic field H. The exponents are defined as follows: 

C = -T?£ = t-°; m = -^ = t^ X=^ = t^, m(T = T c )=H^, (2) 

where t = (1 — T/T c ) is the reduced temperature, C is the specific heat and m denotes the mag- 
netization. These exponents 2 are related (Baker 1990) via, inter alia, the Rushbrooke inequality 
(Rushbrooke 1963) 

a + 2/3 + 7 >2. (3) 

This inequality was originally derived for a second order phase transition. However, it is valid for 
transitions of all orders. For a higher order transition, the singular derivative is the pth derivative 
of the free energy with respect to temperature or magnetic field. By definition, the exponent of 
the thermal derivative is ji. By construction, p — ji = 2 — a and Rushbrooke inequality becomes 
fi + 2/3 + 7 > p. But this is really not yet the relationship we are seeking, i.e., the one between 
the exponents of possibly singular observables. In order to focus on that relationship it is easier to 
study the equality. 

Thus in a pth order phase transition, the singular derivative with respect to magnetic field (as 
an example of a mechanical variable) is also the pth one, N p = d p F/dH p = d p ~ 2 x/9H p ~ 2 . Let 
us consider therefore (while defining Q p and k p ; the exponent k p need not be confused with the 



2 We know of only one book (Pippard 1966) which considers a higher order transition. 
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Landau parameter k above) 

F=-t'-»~-(H C2 -HY>-<>, H C2 =t x \ (4) 

QPp 

N *=8W= ® 
It follows therefore that k p = x 2 ( p . By identifying the most singular terms, we can derive the 
following identities: 

p - n = x 2 (p - (p), (3 = x 2 (p - ( p - 1), S~ 1 =p-C p -1, 7 = -x 2 (p - C P ~ 2 )- (6) 
These can be further processed to produce the following identities: 

P + 7 = x 2 , ( P =P- ^ + + 2 py «p = (P - 1)7 +(P- 2)^9. (7) 

Here x 2 , the exponent for the upper critical field H c2 (T) is an observable, as are (3 and 7. The 
latter two, in general, correspond to non-singular (if p > 2) observables. The first expression 
in Eq. (7) is a scaling law relating temperature dependence of exponents between two different 
experiments. Similarly the last expression too relates different experiments. The middle expression 
is an important algebraic link between the two. 

Finally, by eliminating x 2 we can derive the scaling relations: 

(p - \)n p + pP p + k p = p{p - 1), (8) 

Kp = (3 p [(p - 1)6 P - 1], (9) 

lh=P-Pp(S P + l), (10) 

(8 P + l)Kp + [(p - 1)5 P - l](/ip -p) = 0. (11) 

Recall that the classical values here are n p = k p = 0; (3 P = 1/S P = p — 1 and 7 P = 2 — p. These 
might look unfamiliar to readers more accustomed to critical exponents in magnetism where the 
coupling to magnetic field appears through Zeeman interaction (-M.H). In the case of a gauge 
coupling, appropriate for a superconductor, the classical exponents are exactly as described above 
for p = 2. Equation (8) is the new version of the Rushbrooke (Baker 1990; Rushbrooke 1963) 
equality while Eq. (9) is a generalization of "Widom's law" (Widom 1964). Finally Eq. (10) is the 
well known Griffith's law(Griffiths 1965). Thus, if the exponents are non-classical (p ^ 0) they are 
all interconnected via the scaling laws. 

We have not discussed the exponents v and r] and the associated analog of the Buckingham- 
Gunton inequality (Baker 1990; Buckingham and Gunton 1967,1969). These exponents involve 
spatially varying features and remain subjects for future discussion. However, a calculation of non- 
classical exponents, in principle, is no more difficult than it has been for p = 2. Put differently, 
it is not advisable to consider a free energy with a power law temperature dependent exponent of 
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3.6 as 2 + 1.6, which would describe a second order phase transition with a specific heat exponent 
a = —1.6, due to fluctuations. Indeed much of the machinery for a calculation of exponents is 
designed with small a in mind. More likely, the right course of action is to consider a different 
mean field theory followed by a consideration of fluctuations about that theory, i.e., the free energy 
exponent should be viewed as 4 — 0.4, with the mean field corresponding to the higher integer. 



3. Other Thermodynamic Identities 

The above relations were originally derived to relate exponents which incorporated fluctuation 
effects around a second order transition. In the following we consider other identities which are 
characteristic of a mean field description for a second order transition. These can not clearly 
be extended to include fluctuations but extension to higher order transitions is possible. For 
example near T c , it is possible to derive a relationship between the specific heat (in its temperature 
dependence) and the London penetration depth A(T). That a result such as this should exist is 
clear from the way the magnetic field couples to the system. Consider for example the free energy 
for a second order phase transition in the presence of a magnetic field, 

-a t\^\ 2 + b\^\ 4 + c 



V + — A V 

4>o 



+ -L(VxA) 2 . (12) 



Here t = 1 — T/T c ; a , b and c are all positive constants. The complex order parameter ip is the 
independent variable which takes its ground state value V'o by minimizing the free energy F. A is 
the vector potential so that the magnetic induction B is given by B = V x A. Finally, <fio = h/2e 
is the superconductor flux quantum and (1q is the magnetic permeability. If we minimize Eq. (12) 
with respect to the vector potential to derive the Euler-Lagrange equation for A(r), we get 

h-^£) ]U - (13> 

Thus, near T c , all temperature dependence in A -2 comes from the temperature dependence of the 
order parameter IV'I 2 - The coupling between the vector potential A and the order parameter is the 
gauge coupling and it is clear that this is the only temperature dependence possible for A~ 2 (near 
T c , the quasiparticle contribution to A -2 is expected to be small compared to the order parameter 
contribution here). The free energy F(T) can be written as 

F(T) = -b^ ; V> 2 = ^, *>0 (=0, t<0). (14) 
Thus, by expressing V'o in terms of A, a simple result emerges, 

C(T)- T d2F - bT * \~ 4 (i 5 ) 

Apart from the constant b which appears in C(T) and c, which is a measure of the gradient coupling, 
the other constants such as 4>o and //q are fundamental. 
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Equation (15), to our knowledge has not been widely used. It is specific to the form of the 
gauge coupling and it is easy to generalize to a pth order phase transition. The free energy is 
believed to be: 

2m 



F{T) = -a t\^ p ~ l) +b\^\ 2p + c 



V + —A %b 
00 ' 



(p-i) 



2 



(16) 



where we have omitted the magnetic field energy term, the last term in Eq. (12), for brevity. 
Here the singular derivative of the free energy is M p = d p F/dT p = d p ~ 2 C / dT p ~ 2 , rather than the 
specific heat itself. Thus [I = p/(p — 1)]; 

M b f^Y 9 X~ 21 (17) 

Mp ~ (p-l)(2^c)' {&) dT pX ■ (17) 

Note that for p = 2 Eq. (17) reduces to Eq. (15). 

It is tempting to wonder whether Eq. (17) provides a scaling relation between the exponent 
for free energy (or specific heat) and one for the London penetration depth A(T) (A~ 2 ~ t Xl ). Our 
derivation so far is limited to the mean field regime. Let us define scaling exponents so that near 
% 

F(T) ~ -t p ~^ H C1 = t x \ H C2 ~ t x \ (18) 

then we expect from H 2 = H Cl H c . 2 , x\ + X2 = p — \i while from Eq. (17), we have x\l = p — \i so 
that 

x 1 = ^^(p-l), x 2 = ^. (19) 
p p 

For p = 2, x\ = X2 = 1 — n/2. For p = 4, the classical values are x\ = 3, xi = 1. 

We have measurements (Kumar, Hall and Goodrich 1999; Hall et. al. 2000) for both x\ 
and X2'. x\ = 3.03 ± 0.16 and X2 = 1.21 ± 0.02. However, there are problems. From Eq. (19); 
we see that, in order to be consistent with the temperature dependence of H C 2, X2 < 1. Another 
manifestation would be x\ + X2 < 4; but the measured values do not satisfy that. This is nearly 
but not quite satisfied by the experimentally measured exponents. 

This may well be due to the fluctuation effects. The London penetration depth is a measure of 
the superfluid density. But the temperature dependence of the superfluid density is better defined 
via the penetration depth. When we replace the order parameter in Eq. (14) in favor of A by using 
Eq. (13), we have restricted the identity to mean field theory. Moreover, x\ is determined from the 
temperature dependence of the lower critical field H c \ (T) which contains logarithmic corrections to 
the temperature dependence of A. In view of the divergent temperature dependence (Kumar, Hall 
and Goodrich 1999; Kumar, unpublished 2001) of the Landau parameter n = A/£, the exponent 
for A may be smaller than x\\ £ denotes the coherence length. We have here a relationship between 
the specific heat C(T) and the London penetration depth A(T). An objective of the discussion here 
is to motivate precise and direct measurments of A(T) 
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4. Transition in a Film and Melting of vortices 

A feature of the Ginzburg-Landau model for a pth order transition described by Eq. (16) is 
the temperature dependence of the superfluid density p s oc l^j 2 ^ 1 ) oc t v ~ x . The superfluid density 
may be identified by the kinetic energy being \p s v1 where v s = ^-(V<^>), 4> denotes the phase 
of the order parameter and m* is the effective mass of the Cooper pair. One consequence is the 
nature of the Kosterlitz-Thouless transition (Tinkham 1996) in a thin film with thickness d. In two 
dimensions, a vortex-antivortex pair can unbind above a temperature Tkt-, leading to s destruction 
of the superconducting state. The expression for the transition temperature Tkt is given as 



32,r 2 XHT KT ) 



or 



where d^ is a scaled length. There are two major consequences. Tkt in general is smaller for a 
higher order phase transition than for a second order one. But the ratio Ap s /TxT, a universal 
constant, remains the same constant here. The decrease in Ap s is entirely due to the decrease in 
Tkt- 

In three dimensions, the analog of K-T transition is flux-melting, (Tinkham 1996; Houghton, 
Pelcovits and Sudboe 1989) particularly in anisotropic, high T c superconductors. In the H-T 
phase space, the melting of vortices occurs at a phase boundary which has been calculated by 
among others Nelson (1988) and Houghton, Pelcovits and Sudboe (1989). The qualitative results 
can be encapsulated in a simple expression (Tinkham 1996) for the flux melting line. 

B m (T) ~ A" 4 . (21) 

Given the strong temperature dependence of A, this field cannot be identified with the observed 
irreversibility line. We will defer this subject to a later exploration. 



5. Discussion 

The principal result we have here is the derivation of scaling laws, Eqs. (8)-(ll), appropriate 
for the exponents in a higher order phase transition. In addition, we have explored the possibility 
of other identities which at the moment seem to be restricted to a mean field description but may 
lend themselves to a more general analysis. 

This is in fact related to a more general issue in superconductivity. In the conventional su- 
perconductors, the transition (to the extent known) was second order, and a mean field theory 
was quite sufficient, the exponents were all classical. Then came the high T c superconductors and 
dominance of fluctuations. Thus, in high T c superconductors, it is often argued (Blatter et al 1994) 
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that there is no H C 2(T). The coherence length is very small and the fluctuations are dominant. The 
transition is determined by whether the vortices are pinned in which case the resistivity is zero, 
or not in which case the system becomes a normal metal. This is qualitatively different physics 
where macroscopic defects play a critical role. Scaling laws, developed to describe moderate effects 
of fluctuations, have apparently no relevance in this case. 

But in BKBO the coherence length, as determined by H&iT) is of order 6 nm. This is not small. 
That the transition is higher order does not necessarily mean that fluctuations have altered the 
physical landscape. On the contrary, based on the currently available information, the fluctuations 
here may well be moderate and scaling laws will likely be useful. 

Doubts have been raised with regards to whether the transition in BKBO is of higher order. 
For example Woodfield et. al. (1999) have reported a small discontinuity in the specific heat at the 
superconducting T c . If there is a discontinuity in the specific heat then there is no need to invoke 
a higher order phase transition. But the order of the transition (Kumar, Hall and Goodrich 1999) 
is determined by the temperature dependence of the thermodynamic critical field. The Berkeley 
results (Woodfield et. al. 1999) also indicate that the discontinuity disappears at magnetic 
fields of order 3 Tesla. This is much smaller than any other measurement of upper critical field 
(see other citations in Hall et. al. (2000)). The case developed in Kumar, Hall and Goodrich 
(1999) is internally consistent and is based on measurements of field and temperature dependent 
magnetization. 

Moreover, in a wide ranging study of a number of high T c materials, Junod, Erb and Renner 
(1999) suggest that the specific heat anomaly at the superconducting transition in Bi2212 com- 
pounds is distinctly different from the mean field behavior in conventional superconductors as well 
as a A-point anomaly in Y123 (YBa2Cu307). In particular, they suggest that the specific heat in 
Bi2212 appears similar to that near a Bose-Einstein condensation (BEC). In as much BEC (in an 
ideal Bose gas) may be seen as a third order phase transition 3 , Bi2212 should be considered as 
a candidate material as well. As Junod, Erb and Renner (1999) note, the specific heat in Bi2212 
is continuous. It is the temperature derivative of the specific heat which is singular at T c with an 
exponent of (13 = 0.33. 

Finally, we note that a third order phase transition has been proposed in the large- N limit of the 
two-dimensional U(iV) lattice gauge theory with variation in the coupling constant (or analogously 
the temperature in statistical mechanics) (Gross and Witten 1980). A discussion of the third order 
transition in an associated chiral model and pertinent exponents is given in Campostrini et. al. 
(1995). A possible third order phase transition in Invar type alloys (Shiga and Nakamura 1990) 
and a fourth order transition in the antiferromagnetic Blume-Capel model (Wang and Kimel 1991) 
have been reported in the literature. However, no attempt was made to either provide a free energy 



3 The compressibility of a neutral Bose-Einstein condensate is infinite. Thus London (1964) has argued that the 
BEC in an ideal Bose gas in the P-T plane is a first order phase transition. But in a non-ideal system, we might not 
have this anomalous behavior. 
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or derive the scaling exponents. 

6. Acknowledgment 

This work was supported in part by the U.S. Department of Energy and in part by the National 
Science Foundation. We are grateful to R. Goodrich, D. Hall, D. Hess, A. J. Houghton, G. Stewart 
and J. D. Thompson for discussions. 

REFERENCES 

Baker Jr., G. A., 1990, Quantitative Theory of Critical Phenomena, (Boston: Academic Press). 

Blatter, G., Feigel'man, M. V., Geshkenbein, V. B., Larkin, A. I. and Vinokur, V. M., 1994 Rev. 
Mod. Phys. 66, 1125. 

Buckingham, M. J. and Gunton, J. D., 1967, Phys. Rev. Lett. 20, 143; Also see 1969, Phys. Rev. 
178, 848. 

Campostrini, M.,Rossi, P. and Vicari, E., 1995, Phys. Rev. D 52, 395. 

Ehrenfest, P., 1933, Proc. Amsterdam Acad. 36, 153. This paper is available in 1959, Paul Ehrenfest: 
Collected Scientific Papers edited by M. J. Klein (Amsterdam: North Holland). 

Fisher, D. S., Fisher, M. P. A. and Huse, D. A., 1991, Phys. Rev. B 43, 130. 

Griffiths, R. B., 1965, Phys. Rev. Lett. 14, 623. 

Gross, D. J. and Witten, E., 1980, Phys. Rev. D 21, 446. 

Hall, D., Goodrich, R. G., Grenier, C. G., Kumar, P., Chaparala, M. and Norton, M. L. 2000, Phil. 
Mag. B 80, 61. 

Houghton, A. J., Pelcovits, R. A. and Sudboe, A., 1989, Phys. Rev. B 40, 6763. 
Junod, A., Erb, A. and Renner, C. 1999, Physica C 317-318, 333. 
Kumar, P., Hall, D. and Goodrich, R. G. 1999, Phys. Rev. Lett. 82, 4532. 

Kumar, P., "Theory of a Fourth Order Phase Transition: Superconducting Transition in BKBO", 
UF preprint (2001). 

London, F., 1964, Superfluids, V II: Microscopic Theory of Superfluid Helium, (New York: Dover), 
p.52 



-10- 



Methfessel C. E., Stewart, G. R., Matthias, B. T. and Patel, C. K. N., 1980, Proc. Natl. Acad. Sci. 
USA 77, 6307. 

Nelson, D. R., 1988, Phys. Rev. Lett. 60, 1973. 

Pippard, A. B., 1966, Classical Thermodynamics, (Cambridge:Cambridge U. Press). See p. 142 and 
the discussion following p. 146. 

Rushbrooke, G. S., 1963, J. Chem. Phys. 39, 842. 

Shiga, M. and Nakamura, Y., 1990, J. Mag. Mag. Mat. 90-91, 733. 

Tinkham, M., 1996, Introduction to Superconductivity, Second Edition (New York: McGraw-Hill). 
Wang, Y. L. and Kimel, J. D., 1991, J. Appl. Phys. 69, 6176. 
Widom, B., 1964, J. Chem. Phys. 41, 1633. 

Woodfield, B. F., Wright, D. A., Fisher, R. A., Phillips, N. E. and Tang, H. Y., 1999, Phys. Rev. 
Lett. 83, 4622. There was an earlier report as well, Graebner, J. E., Schneemeyer, L. F. and 
Thomas, J. K., 1989, Phys. Rev. B 39, 9682. The anomalies are small and are obtained after 
considerable processing of the data. 



This preprint was prepared with the AAS IATjrjX macros v5.0. 



